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ABSTRACT 

This document applies the cylindrical waves associated with the Sommerfeld solution for 
a diffracting half-plane in the evaluation of the effects of secondary radiation on the focal- 
region-fields of paraboloids. A cross-correlation between the computed focal-region-fields and 
the feed-antenna-fields is then introduced to determine the modified radiation pattern of para- 
boloidal antennas. The half-plane solution is adapted to estimate the effects of a polygonal 
fence on 9-meter S-band and 3-meter S-and Ku-band antennas at the Merritt island, FL Space- 
craft Tracking and Data Network Station. 
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A RADIO-FREQUENCY ANALYSIS OF PARABOLOIDAL ANTENNAS 
LOCATED NEAR DIFFRACTING FENCES 


INTRODUCTION 

This document is an attempt at estimating the effects of proximate diffracting fences on 
the performance of paraboloidal antennas at tiie Merritt Island, FL Spacecraft Tracking and 
Data Network Station. Both 9-meter S-band and 3-meter S-and Ku-band installations are of 
concern. Mathematical simulation of the paraboloids and the diffracting fences affords a means 
of assessing the effects as a function of antenna and fence parameters. High-speed digital com- 
puters may be used to extract quantitative information from an analysis based on the vector 
Kirchhoff-Kottler physical optics theory for continuous surfaces and the Sommerfeld diffraction 
theory for half-planes. 

A plane-wave which is incident on a paraboloidal antenna induces an electromagnetic 
sheet- current on the latter. From this a highly satisfactory computation of the focal-region 
fields, wavefronts, and time-average Poynting vectors can be made. The assumed primary wave 
incident on the paraboloid .may also illuminate diffracting fences which act as secondary 
sources. Each fence will radiate a cylindrical wave exhibiting, in general, a linear phase grad- 
ient, and a non -trivial field-intensity distribution. If phase is preserved in a mathematical sim- 
ulation, and all polarization “components are included, the total perturbation of the sheet- 
current on the paraboloid may be calculated and, subsequently, a new set of focal-region fields, 
wavefronts, and time-average Poynting vectors may be determined. 

The far-field radiation pattern of a paraboloid, influenced by one or more diffraction 
fences, may in principal be found from a knowledge of the focal-region fields and feed-antenna 
fields over some surface of convenience. A sphere centered on the focal-point of the para- 
boloid, with radius equal to the far-field distance of a feed horn, is satisfactory for this pur- 
pose. A cross-correlation of the focal-region and feed antenna fields, by means of an algorithm 


that is responsive to the spatial (vector) and temporal (complex-phase) character of the prob- 
lem leads to a single voltage for every angle of arrival of an assumed plane-wave, In this man- 
ner, the far-field or transmission radiation pattern evolves, from an analysis which is predicated 
entirely on reception arguments, in the presence of diffracting fences, 
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PLANE-WAVE SOURCES 

Two orthogonal plane-wave sources are admitted in the simulation. See Ref. 1, p. 578, 

E - polarization 

E -• (- cosa sini3» - sinasln^, cos|3) (1) 

H “ (- sina, cosa, 0) e"*’'® (2) 

H - polai/ization 

E = (sina, -cosa, 0) e"*^^ (3) 

H != (-cosasin/3, -sinasin/3, cos0) e"**'^ (4) 

The phase factor of the plane waves is 

g-iks = g-ik (x cosa cos/3 + y sina cosj3 + z sin/3) _ (> 5 ) 

Angles (a, (3) and coordinates (x, y, z), above, were originally in a local context for a 
half-plane analysis as presented by Born and Wolf, but are now utilized in a global context 
(GSFC coordinates) for the incident plane wave of the simulation. A given point (x, y, z) 
on the paraboloid will appear as (x", y", z") at the origin of a rotated and translated diffrac- 
tion fence. Similarly, angles (a", /3") will emerge as the local angles of arrival at such a fence. 

See Fig. 1. It can be seen that (a) is an azimuthal variation in the equatorial (x, y) plane, 
and (P) is a polar variation emanating from the (x, y) plane, 

Any degree of polarization ellipticity, orientation, and rotation sense of the polarized state 
is available via either (1) and (3) or (2) and (4). It may be more direct to combine two 
counter-rotating circular states, in any given instance, since the relative phase of the circular 
states and inclination angle of the polarization -ellipse are simply related, 

= 21//J , (6) 
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the axial ratio of the ellipse derives from the relative weights, 


AR = (wj + co2)/(coj - CO 2 ) 


( 7 ) 


and the rotation sense of the field vector is determined by the la:;;er of (^i) and (CO2). 


OttlaiNAL 
OF POOR 

SCATTERING FORMULATION (PARABOLOID) 

An amended form of the compIex*vector Kirchhoff-Kottler formulation (physical optics) 
suffices to compute the focal region fields, wavefronts, and time- average Poynting vectors of 
the paraboloid. It is a variation of the Stratton-Chu formula, See Ref. 2, p. 460; Ref, 3, 
p, 158; Ref. 4, p. 4; Ref. 5, p. 35. 


PAQE 

QUAUTY 


R(x, y, z) = 


1 1 


jcoe 47 t f j (n x Hj) , (8) 

Sji Sj 


H(x, y, z) “ ” / (n X Hj) X 

^|; = / r , 


ds 


(9) 

( 10 ) 


where 


mh 


- + (jk = — ) T. 

r ‘ 


(ID 


The phase of the Green's function (10) is seen to be consistent with the plane-wave formula- 
tion of Born and Wolf, given previously. 
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SCATTERING FORMULATION (HALF-PLANE) 


ORIGINAL PASEJI 
OF POOR QUALfTY 


The total fields, in local coordinates, for diree-dimensiomil diffraction of a plane wave 
by a half-plane are given in Ref. I , p. 579, in the Gaussian system of units as 
(E-plane polarization) 

- -Hy sin 13 (12) 

Ey « sin 13 (13) 

E^ sj c"7ri/4 gos d e**^^*^ cos - z sin /?) [G(p) - G(q)] (14) 


Hjj K - e ^ ^ < sin£v:[G(p) + G(q)] -H 1 y 

7 ^ ( 


/ 2 


kr cos p 


sinoi/2 cosd/2| 
(15) 


K 


0 “^ i/4 Q ik(r cos /3 - zsin P) |ccs«[G(p) - G(q)j - ^ V sina/2 cosO/2 1 

— ( i cos p ) 


(16) 


H, H. 0 


(17) 


As mentioned previously, (x, y, z) above will be redesignated (x", y", z") due to certain 
correspondences, translations, and rotations among tlie Cartesian coordinate systems used hereir 
The companion expressions to (12) through (17) for the H-plane polarization case are not 
given explicitly in Ref. 1, It happens that the E-plane polarization case was developed by 
means of two auxilliary equations, 


and 


/ - 1 sm ^ 9U - 1 sm 9U 
V k 9x ’ k 9y ^ 


(18) 




iJ i. JH 

\ k 9y ’ k 9x 


, 0 


. (19) 
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The H-plane polarization ease may, therefore, be developed using two other auxilliary 
equations. See Ref. 1. p. 579. 

i 3U - i 3U 


^ = (i 


k 3y ’ k 9x 


. 0 


( 20 ) 


and 


H " 


i sin /? 9U - i sin j3 9U 


9x ’ 


9y 


, cos^ /3 U 


( 21 ) 


Details of the derivation may be found in Appendix A of this document. The results are 
given below, also in the Gaussian system of units, as 
(H-plane polarization) 




e“7d/4 gik(r cos jS - zsin /3) | 

7 ^ 


f 2 a 0 

sino![G(p) - G(q)] © iV- — cos sin 

^ kr cos p 2 2 


( 22 ) 


E., = 


Q- 7Ti/4 gik(r cos /3 - zsin ^) |coso;(G(p) + G(q)l © i\/^- — cos cos — | 

— = I kr cos )3 2 2' 

VT 

(23) 

Ej, = 0 (24) 


= Ey sin |3 

Hy = - Ej^ sin p 
_ £^i/4 glk(r cos p - zsin P) [Q(p) 
" " Vlr 


(25) 

(26) 
(27) 


Note; 

Present derivation (+ i) in equation (22), (+ i) in equation (23). 
Ref. 1, 2nd Ed. (- i) in equation (22), (- i) in equation (23). 
Ref, 1, 6th Ed. (- i) in equation (22), (+ i) in equation (23). 
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The evaluation of G(p) and G(q) presents alternatives. Assume that 

/J < 7T/2 


(24) 

(25) 


See Fig. 2, where the diffracting half-plane is bounded by the ±z axis and lies in the (+x) 
part of the xz-plane. It can be seen that 


and 


P > 0 I 


q > 0 I 


6 > IT + a 


0 > 7T - a 


(26) 


(27) 


from an inspection of equations (2A). That is, p > 0 in the shadow region (SR) and q < 0 
in the reflection region (RR). 

It can also be seen that the condition for asymptotic evaluation, 


kr » 1 , (28) 

given in Ref. 1, p. 572 for the two-dimensional diffraction of a plane wave by a half-plane, 
becomes 

kr cos |3 » 1 (29) 

for the present three-dimensional case. Further, (p) and (q) may vanish under same conditions: 


and 


P 


0 


0 = rr + a 


(30) 


?.t the shadow and reflection boundaries, respectively, even for large (kr cos /3). This indicates 
that the requirement on (a) for asymptotic evaluation implies 


p » 1 , q » 1 
8 


( 31 ) 


From Ref. 1, p. 569, 
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G(a) = F(a) (30) 

and parameter (a) may be zero, positive, or negative. Here, 

F(a) = / e‘/^^d/Li . (31) 

a 

The development of Ref. 1, p. 572 pertaining to the evaluation of G(a) for negative (a) is not 
explicit for small values of (a). Appendix B of this document gives the details of a develop- 
ment leading to a general formula. 


G(a) = e"^^^ [ yJlF U(- a) + sgn(a) / e*/^ d/a ] 

13 | 


(12), (7-B) 


where 


and 


U(x) = 1 1 X > 0 
0 i X < 0 

sgn (x) = + 1 j X > 0 

- 1 1 X < 0 


(33) 


(34) 


For large values of (a), equation (32) above reduces to 

G(a) = y/lT U(- a) + sgn(a) + 0 — , (35) 

2|al ia|3 

2 

which agrees with Ref. 1, excepting a typographical error; (e"^^) should read (e"^“ ), equation 
(32), page 572. 

See Appendix C for a conversion of the Ref. 1 equations, above, to the rationalized MKS 
system of units. 
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COORDINATE SYSTEMS 

A convenient means for computing the sheet-current induced on the paraboloidal reflector 
by the diffracting half-planes or fences is to retain the local coordinate conventions of Born 
and Wolf, Ref, 1, pp. 556-580. Figure 3 illustrates a reference coordinate system (x, y, z) of 
the overall problem and the local coordinate system (x', y', z') of p half-plane shown translated 
from 0 to O'. It can be seen that a correspondence exists between the two systems prior to 
translation of origin. 

(X, y. z) (-x', z', y') . (36) 

The Cartesian components of vector (r') from 0' to a point on the paraboloidal reflector 
( 7 ) may be specified in the (x, y, z) coordinate system using 

r' = r - rjp = (x, y, z) (37) 

and rewritten as 

r' = (x', y', z') (38) 


using the correspondence, (36). 

Prior to rotation of the half-plane, the Born and Wolf local cylindrical coordinates 
(r', 0 ', z') are then found from 


or 


and 


x' = r' COS0' (39) 

y' = r' !?in0' (40) 

(r')2 I r- |2 = (x')2 + (y')2 . (41) 
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Rotation of the liaif-plane results in a re-designation of the local coordinate frame as (x", 
y", z") with origin (0") at (O') as shown in Fig, 4. The Cartesian components of vector (f') 
are now different from those of (r"). The latter are found using 

(x", y", z") = [A] Cx', y', z') (42) 


where 


[A] 


(cos7cosa - cos/3sino:sin7) (0QS7sinQ: + cos/?cosasin7) (siii7 sin/3) 
(-sin7coso: - cos/3sinacos7) (-sin7sina + cos/3cosacos7) (cos7 sin/3) 
(sin/3 sina) (-sin/3 cosa) (cos/3) 


(43) 


See Ref. 6, p. 107 for a discussion of Eulerian angles. 

The frame (x", y", z") is visualized as rotating with respect to the frame (x', y', z'), 
and requires the rewriting of (39) (40) and (41) in double-prime notation to obtain the Born 
and Wolf local cylindrical coordinates (r", 6", z"), Once the latter are obtained, it is conven- 
ient to drop the double-prime notation and enter the Born and Wolf equations with the 


derived values, 


POLYGONAL FENCES 

Additional complexity arises when the half-plane formulation is applied repeatedly to 
mathematically simulate a polygonal fence such as the one illustrated in Fig. 5. Also see flow 
charts, Figs, 6 and 7. The procedure is justified by the parameters of the present problem 
and the inherent behavior of the half-plane formulation when kr » 1, Only S- and Ku-band 
frequencies are to be introduced into the present calculations, and the fence height and length 
for each polygonal element exceed 8-feet. Under these conditions a segment of a cylindrical 
wave is associated with each segment of the fence. Ref. 1, p. 572. Also see Fig. 1 of this 
document. Points (?', r'') which either exceed the cylindrical bounds of a finite cylindrical 
wave, or imply intersection with paraboloid (7^), are excluded from consideration. 

Assuming that a single plane-wave illuminates the paraboloid directly, together with the 
several elements of a polygonal fence, it can be seen that a distinct set of input parameters 
must be associated with each element of the fence. The local plane wave angles-of-arrival 
(a, ^), the local Eulerian angles (a, j3, 7), and the displacement vectors (fj^p) are required. 

A direct method for finding the set of plane-wave angles-of-arrival is to proceed on the 
assumption that the fence geometry, in relation to the paraboloid, has been established. Sel- 
ection of a single pair of values (a, /3), for purposes of illuminating the paraboloid, then fixes 
all of the local angles of arrival for the half-planes. If a unit vector with Cartesian components 
in the (x, y, z) reference frame is chosen co-linear, but antidirectional with respect to propaga- 
tion vector of the plane-wave so that 

Rj = (x, y, z) = (cosa cos/3, sina cos/3, sinjS) , (44) 

the wave angles of arrival are obtainable by combining the correspondence (36) and the rota- 
tion (42). The translation (37) may be ignored here since only angles are of concern. In each 
instance the Eulerian angles of the rotation matrix are known and a vector (x", y", z”) is the 
output. The local (a, /3) for each half-plane is then found by solving a system of equations 
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similar to (44) of the present document. Note that R| is different between Fig. 2 and Fig. 3 
due to the preservation of local (a, P) and [A] not equal to the identity matrix. 

Assume that selection of the (a, P) pair taken with respect to the (x, y, z) frame leads to an 
incident linearly-polarized magnetic field whose vector components, without regard to phase, are 

H = (H^, Hy, H^) . (45) 

via either equation (2) or equation (4), This field must now be decomposed into an EpQL 
and HpQL field, so that, in general, the magnetic fields scattered by each half-plane may be 
evaluated via equations (19) and (21) for E-plane and H-plane polarizations, respectively. It 
follows, using (2) and (4), that the field of (45) may be written in local coordinates as 

(Hj^", Hy", H^") = mj(-sina, coso:, 0) + n|(~cosa: sin/?, - sina sin/3, cos/3) , (46) 

where (mj) and (n^) are the unknown weights of two orthogonal, cophased plane-waves incident 
on the half-plane. The angles (a, /3) of (46) are now the local angles of arrival, computed 
previously, so that (46) is a system of three equations in two unknowns. This leads to 

(Hjj", Hy", H^,") = nij H£_pql + nj Hj^_pQL . (47) 

In the event that an elliptically-polarized phane wave is incident on the paraboloid ( 7 ^) 
and the half-planes or fences, a second and similar resolution of the type discussed is required. 
This leads to an illumination of the form 

HpOT - (^1 Hfi-POL ”1 %-POL^ '^^2 %POL "2 ^H-POL^®^'^2 

where (mj, nj, m 2 , 02 ) are scalar weighting factors and 1 // 2 ) are scalar phase-constants 
affecting rotation sense, axial ratio, and inclination angle of a preferred elliptical state. 

A correction remains to be made concerning phase. The solutions of the half-plane prob- 
lem were originally obtained using a phase factor given by (5), with the origin of coordinates 
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on the diffracting edge. Since the incident plane-wave has been referenced to the origin (0) of 
the (x, y, z) frame the phase at the various origins ( 0 ") is given by a factor 

F = (49) 

which now must be applied to the scattered fields of the half-planes or fences. The are, 
in general, distinct when several fences are encountered, but there is only one value of Rj 
(antidirectional to the propagation direction of the incident plane wave). See Fig. 5 which de- 
picts a polygonal fence. 

The total fields associated with each half-plane or fence are evaluated in a local context" for every 
point on the paraboloidal reflector ( 7 ^) and must be converted to the global context of the problem. 
Since the magnetic field components are expressed in the (x", y", z") system of coordinates, 

(H^', Hy', H,') = [A]“i (H^", Hy", H^") . (50) 

Once (Hx^ Hy^ h^ve been determined, the correspondence of (36) is used in the inverse sense. 

(H^',Hy', H^') (-H^.H,, Hy) (51) 

The translations (F|^) are irrelevant since translation does not affect projections between two 
non-congruent coordinate frames. In the preceding [A]"^ is, in general, distinct for the half- 
planes or fences since the Eulerian angles (a, /3, 7 ) are distinct. 

Once the total magnetic fields associated with the diffracting fences have been computed and 
the geometric optics field has been subtracted for each fence, a superposition of sheet currents 
may be visualized. That is, 

Ktot “ 2h x H. + 2n X (Hsj + HS 2 + Hs^ + ...Hsj) (52) 

(primary induction) (secondary induction) 

for (i) half-planes or fences. The focal region fields may then be computed using equations ( 8 ) 
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through (11) of this document in conjunction with the incident magnetic fields of (52). See 
Figs. 6 and 7 for a flow chart of the preceding operations, and Appendix D for a discussion 
relating to the total field. 
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CROSS-CORRELATION (FOCAL-REGION FIELDS AND FEED FUNCTION) 

An algorithm is now written to express the received voltage as a cross-correlation between 
a completely polarized focal-region and a feed function. Such a voltage represents one point of 
an aperture-antenna radiation pattern, The cross-correlation is effected on a spherical locus, 
as shown in Fig, 8, for convenience. Its radius is the far-field distance of the feed horn, 

Rs = , (53) 

A 

The algorithm derives directly from the fact that the induced voltage which results from 
the interaction of two elliptical states is given by 

Vqc = d* • P . (54) 

where q* is the complex-conjugate of the received field and ^ represents a linear receiving aji- 
tenna. See Ref, 8, p. 105, Ref. 9, p. 184. When p and q are both unit complex-vectors, 
an expression for polarization efficiency may be written as 

hpoL ^ I P • q'*' P • (55) 

In the subject cross-correlation for receiving apertures equation (54) becomes 

Vqc = * E/ dl2 . (56) 

Here (E^^ is a complex-vector focal-region field and (E^) is characterized by the feed function. 

In the present context the received field may be regarded as having the dimensions (volts/meter^ ), 
and the feed function acts as an effective height (meter^) per unit solid angle (S2). The dif- 
ferential solid angle (dfZ) may be formally associated with area (dS) on the spherical locus (S) via 

dn = dS (57) 
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IF I s Rs 


( 58 ) 


and the unit normal (iTj) is that of a sphere. See Ref. 10, p, 290. It is noted that (Ej^ and 
(Ey>) are, in general, not unit quantities in (56) and only relative voltage is required for 
a radiation pattern. 

The magnetic counterpart to (56) may therefore be written as a proportionality: 


Vqc a // Hr* • H; cia . (59) 

The magnetic field (Hj^ incident on tiie locus (S) is obtained from equation (9). The term 
(H/) depends on the feed function selected for the calculation, It may be as simple as 



A cos^ CB) 
P 




(60) 


where 


hi 


p X V 
I p X V 1 


(61) 


Here (A) is a scalar amplitude, (N) is a directivity factor, CH)is a polar angle taken with respec 
to the feed horn axis, (v) is a polarization moment, and 


IpI = Rs . (62) 

Many other representations are possible: rigorously derived expressions for feed horns and open 
waveguides, numerical data, etc. 

Returning to Fig. 8(a), it appears that a rotationally symmetric situation exists when the 
primary or incident plane-wave is axially directed (a = tt/ 2, /3 = 7r/2). This is somewhat illusory 
since the rule for current induction 

K = 2n X Hj (63) 
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converts the axial symmetry associated with the paraboloid to uniplanar symmetry, This uni* 
planar symmetry persists for (a < tt/2, « tt/ 2), or paraxial wave arrival, although the received 

fields (Ej^ are now displaced with respect to the feed function (E^), as suggested by Fig. 8(b), 
Introduction of diffraction fences, Fig. 8(c), will leave the calculation entirely without sym- 
metry, 
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CONCLUSION 

This document constitutes tiic first or analytical phase of an attempt to estimate the ef- 
fects of proximate diffracting fences on the performance of paraboloidal antennas at the Mer- 
ritt Island, FL Spacecraft Tracking and Data Network Station, The second and third phases 
are the writing and verification of a computer program which incorporates the theory presented 
herein, and tlie actual generation of numerical data for assumed antenna and geometrical para- 
meters. Introduction of fence mesh details in a higher level analysis, attendant program modifi- 
cation, and subsequent numerical calculations, constitute possible additional phases, 

An appreciation of the antenna environment may be formed from Fig, 9, which shows 
the Kennedy Space Center Layout (KSC). The effects of the diffraction fences are seen to be 
a function of antenna pointing angle, elevation, range, and frequency. A photograph showing 
the 30' (9-meter) USB Antenna, No, 2, is included as Fig. 10 and should be compared to the 
abstraction originally presented as Fig. 1. It appears, from the large number of parameters, 
the repeated appearance of integrals in the preseiit formulations, tlie relatively high frequencies 
of operation (S- and Ku-bands) and the electrical size of the antennas (3-meter, 9-meter), that 
a highly organized approach to the overall problem is Imperative. Accordingly the programming 
is proceeding with careful monitoring of computer core and cpy requirements. Option switches 
are being incorporated as suggested by Figs. 6 and 7. Initial computations will be restricted 
to linear polarization and single-fence configurations to obtain program verification although 
the program will be written to accommodate arbitrary angle of primary wave arrival, elliptical 
polarization states for the received fields and the feed horn characteristics, and polygonal fences 
with elevation contours included. 
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Figure 1. Source (PW) and Secondary 
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Figure 2. The Shadow and Reflection Regions (SR, RR) for a > rt/2, /3 
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Figure 6. Flow-Chart Showing Computation of Focal-Region Fields (FRF I) 
and Perturbation (FRF II). 
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Figure 7. Flow-Chart Showing Computation of Combined Focal-Region Fields 
(FRF III) and Radiation Pattern 
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(b) Paraxial Plane-Wave (PPW) 



(c) Paraxial Plane-Wave (PPW) 
and (Fj) Present 




Figure 8. Cross-Correlation of Feed-Pattern 'E^) and Received-Field (E^) on Spherical 
Locus (S) 
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Figure 9. Kennedy Space Center Layout (KSO 
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^lgure 10 30-Foot (9-M) USB Antenna No. 
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Scattered Fields for H-plane Polarization 

I 

The evaluation of equations (20) and (21) requires the partial derivatives (11,^) and (Uy) 
for the H-plane polarization case, An expression for (U), given in Ref, 1, p. 580, may be used 
to form (Ujj) and (Uy) for the E-plane polarization case only, and these partial derivatives may 
be verified by utilization in equations (18) and (19) to obtain the explicit field equations (12) 
through (17). This process is now carried out in detail. Subsequently an expression for (U) 
of the H-plane polarization case is found and the process is repeated so that explicit field 
equations may be written. 

E-piane Polarization: 

n = sec/Je^^^C*^ -os/J - zsin/3) fQ(p) „ = ce^krcos/3 [qq,) _ Q(q)] ,= ce^E _ Q^qj] 

s/W 

Here C = sec/J/>/i7 and E = ik cos/3. 

1 I 

p = ~(2kr cosj3)*'^ cos (0 - oi) , q - -(2 kr cos/3)''^* cos (0 + a) 

x = r COS0 , y = r sin0 

_ip2 s cos/J(r + xcosa + ysina) = - i(E + A) 
where A - kcos/3(x cosa + y sina) 


(1-A) 

(2-A) 

(3-A) 

(4-A) 


-iq^ = -ik cos/3(r + xcosa - ysina) = - i(E + B) 
where B = kcos/J(x cosa - y sina) 


(5-A) 


G(p) = e'4P^ / e^^^ d/u 
P 


* 2 1 1 *2 

, G(q) = e“*^ / d/i 


(6-A) 


U = Ce“^^ / e^^ d/i - Ce"’® / e^^^ dp 

P q 

TT ^ "lA d oo . 2 , ^ _iR 9 oo , 2 , ^ 9e 

U„ = Ce T— f dp - Ce - / e‘P dp + C t-- 

^ 8x p 3x q 9x 


A-1 


-lA 


1..2 uc — i,,ti 

f e*P dp - C — f e*P dp 
p 3x q 

(8-A) 


(7-A) 


9e"l®oo 


Applying Leibnitz’s Rule to equation (8-A), 
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“ -Ce"^^ ei(E+A) ^ + Ce"^^ -Cik cos/3 cosa ^e"*^ d/4 ~ e"*® e*^^ 


After simplification, 


U. - -Ce‘E 


ik cos/3 cosa [G(p) - G(q)] + ~ (p - q) 

9x 


where 


and 


p - q = -4=^(2 k cospf'^^ sin •— t (x^ + 

Y 2 2 


~ = iLif gj„ JL sin 

dx " 2 2 


Substituting (12-A) into (10-A), and using 


Hy = 


U. 


equation (16) is verified. 

A similar process produces 


(9-A) 


(10-A) 


(11 -A) 


(12-A) 


(19) 


Uv = ~Ce^^ ! ik cos/3 sina [G(p) + G(q)] + (p ~ q) 


and 


9(P - q) 


3y 


cos 0/2 \ 9(p - q) 


9y \sin 0/2 

Substituting (14-A) into (13-A), and using 


9x 


(13-A) 


. (14-A) 


Hx = -r Uy 


, (19) 


equation (15) is verified. 

An expression for (U), appropriate for the H-plane polarization case may be forced 
directly from a specialized case (/^O) in Ref, 1, p. 574. Using (i'll) of the present document, 
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e-)Tl/4 ik(r cos|3 - zsln|J) Jk r cos/J 

U a _ sec/3 e (G(p) + G(q)] « Ce [G(p) + G(q)] * [G(p) + G(q)] (15-A) 

Using equations (2-A) through (6-A), 


U = Ce”^^ / dfi + Ce”^^ / 6(i 

P q 


(16-A) 


U = Ce"^^ dpi + Ce“*® d/i + C / e^^^ dp. + 


9 « iii 2 


« tii 2 


0x p 


9x q 


C / e*^^ dp 

9x q 


(17-A) 


Using Leibnitz’s Rule, 


U = -Cte-‘A e'(E+A) i£. . Ce-IB e'(B+B) ii 
^ 9x 9x 


- Cik cos/3 cosce ( e“^^ / dp + / e^^^ dp 

V P q , 


After simplification, 


Ujj = - Ce^*® I ik cos/3 cosa [G(p) + G(q)] + (p + q) | 


(18-A) 


(19-A) 


where 


p + q = — ^(2k cos/3)'^* cos — [(x^ + 
v 2 2 


(20-A) 


9 ^ (2k cos/3)'^^ a 6 

— (p+q) = -^ COS- COS- 

Substituting (21 -A) into (19-A), and using 

i 9U 
“ k, 9x 


(21-A) 
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equation (23) is obtained. 

A similar process produces 


3 

Uy “ - Ce^^ ik cos/3 sina[G(p) - G(q)] + (p + q) 

^ 9y 


9 , . (2 k cos/3)’'^ 0! , 0 

— (p + q) = -^ oosysm- 


Substituting (23-A) into (22-A), and using 


equation (22) is obtained. 


L -IH. 

k 0y 


(22-A) 


(23-A) 


Returning to equations (20) and (21), (24) through (27) are obtained immediately using 
(U), (Ujj), and (Uy). Equations (22), (23), and (27) should be compared to those of the spec- 
ialized case (|3=0) of Ref, I, p, 574, 
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Tile Fresnel Integrals 
a > 0; 

F(a) - 

a < 0; 

F(a) = 

Let /i - /u so that 


Wi2 

/ e*^ dfx 

a 

00 . 2 
f e*^ dju 

-|a| ^ 
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jj,/y ■"*0 n 

f dfi = 

^2 = 0 i.,2 

/ e^l^ d/x = 

1^1 , 2 
/ e^l^ d/z 

Ml = -Ia| 

Ml = lal 

0 

A 

o 

°° s.,2 ®® 

e^^^ dfx 


using (2-B), (3-B) and 

e^^^ dju = e‘l^^ d// - e*^^ d/x 

From Ref. 1, p, 573, equation (35), 

G(C) = dM = -~ v/F e^i/4 

so that 

G(a) = [yF U(- a) + sgn(a) /,^| e*^^ dju] 
for 

U(x) = 1 I X > 0 
0 I X < 0 
anh 


(1-B) 

(2-B) 

(3-B) 

(4-B) 

(5-B) 

(6-B) 

(7-B) 

(8-B) 
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5gn(x) - + 1 I X > 0 
- 1 I X < 0 

is a general formula for (a) positive, negative, or zero. 

For G(a) to be continuous at a = 0, 

lim G(- e) = lim G(e) = G(0) = ^ e^‘/4 
e 0 e 0 2 See Ref. 1, p. 76. 


G(- e ) - ^ [s/f - / e‘^^ d/a] 

|e| 


G(e) 


= 0-1 e 


2 oo 


/ e^^ d/x 

Ie| 


lim , « d/u = 

|e|->0 |el 2 


therefore (lO-B) is satisfied. 
Using 


. ^ OO .1 

e~*® / e‘^ d/u 

|a| 


■ + o' ' 


2|a| 


lap 


when (a) is large, 

G(a) = s/W e"‘“^ U(- a) + sgn (a) — y— + 

2|a( \|a|3y 

See Ref. 1, p. 572, equations (29) and (30) for the origin of equation (14-B). 
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(9-B) 

(10-B) 

(11-B) 

(12-B) 

(13-B) 

(14-B) 

(15-B) 
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Units and Dimensions 

The Kirchhoff-Kottler (physical-optics) formulation of equations (8) through (11) of this 
document is in the rationalized MKS system of units. It is characteristic of this system of 
units that, in the far-field, the ratio of the magnitude of the electric to the magnetic field 
yields the free-space impedance (Z^), This is easily seen from 


iE l(v/m) 

I H I (AT/m) 


l-jwMo 47 t 4, I 

I — ]— f (ri X Hi ) X Vi// dS I 
47T Sj ^ 





= Zq (ohms) 


(1-C) 


using (8) and (9) as r ->• oo . 

Tlie plane wave formulation, (1) through (4) and the solution for the half-plane, (12) 
through (17) and (22) through (27), as obtained from Ref. 1 is in Gaussian units. Conver- 
sion of units may be formally effected using Ref. 11. In view of the fact that the magnetic 
fields are of prime interest throughout this paper it is convenient to multiply all electric fields 
of Ref. 1 by a numerical factor 





> 


(2-C) 


Suppressing the units ohms for (Fj^). The electric field and magnetic field expressions obtained 
from Ref. 1 may then be taken as volts/meter (V/m) and ampere-turns (AT/m), respectively. 

In the far-field 


^ ^FF * ^ p * ^V/m * 

I Hff I ” I HaT/f, I 


= Zq ohms 


(3-C) 
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Tlie Total Field 

Equations (12) through (17) and (22) through (27) represent the total fields associated 
with a given half-plane. This may be seen by considering the special case (j3 = 0) and the 
electric field in the reflection region (p < 0, q<0), the illuminated region (p<0, q>0), and 
the shadow region (p>0, q>0), respectively, for a remote observer (r “>oo), The G(p) and 
G(q) of the field equations, cited above, then lead to the indentification of interfering inci- 
dent and reflected plane waves, a transmitted or incident plane wave, and the degenerate or 
zero plane wa'fe for the regions named upon application of half-angle identity formulas. Col- 
lectively these constitute the geometric optics field for the special case being considered: 

t g-ikr cos (6 - a) _ g-ikrcos (9 + a) ^ Q 4 9 < it ~ a 
g-ikr cos (0 - a) , tt - a < 0 < 7T + a 

0 ,7r + o:<0.^27r 

See Ref. 1, p. 572, eqn. (33). 

From the preceding it follows that the geometric optics field is contained in the field 
solution of each half-plane 0 in general), for electric and magnetic fields, and must be 
subtracted (i) times when (i) fences are considered in the present simulation. The paraboloid 
is illuminated directly by a specified incident plane-wave, and will not lie in either the shadow 
or reflection regions. 
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